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The cross section for a Majorana Dark Matter particle annihilating into light 
fermions is helicity suppressed. We show that, if the Dark Matter is the neutral 
Majorana component of a multiplet which is charged under the electroweak in- 
^\ teractions of the Standard Model, the emission of gauge bosons from the initial 

state lifts the suppression and allows an s-wave annihilation. The resulting energy 
spectra of stable Standard Model particles are importantly affected. This has an 
impact on indirect searches for Dark Matter. 



1 Introduction and setup 

The fluxes of stable Standard Model particles that originate from the annihilation (or decay) of 
Dark Matter (DM) in the galactic halo are the primary observable for DM indirect searches. 
The radiation of ElectroWeak (EW) gauge bosons from the final state of the annihilation 
process turns out to have a great influence on the energy spectra of stable particles and hence 
on the predictions for fluxes to be measued at Earth [1,2] (see also Refs. [3,4] for related 
analyses). In particular, there are three situations where the effect of including the EW 
corrections is especially important: 

1. when the low-energy regions of the spectra, which are largely populated by the decay 
products of the emitted gauge bosons, are the ones contributing the most to the observed 
fluxes of stable particles; 

2. when some particle species are absent if EW corrections are not included, e.g. antipro- 
tons from W/Z decays in an otherwise purely leptonic channel; 

3. when the 2 — y 3 annihilation cross section, with soft gauge boson emission, is comparable 
or even dominant with respect to the 2 — y 2 cross section. 

A possible realization of the latter condition has been studied in Ref. [2], where we considered 
the DM as a gauge-singlet Majorana particle \ of mass M x annihilating into light fermions / 
of mass nif <C M x ; it is well known that in this case the XX ~~ ** f f cross section is suppressed. 
In fact, one can perform the usual expansion of the cross section 

va = a + bv 2 + 0(v 4 ) , (1.1) 

where v ~ 1CT 3 is the relative velocity (in units of c) of the DM particles in our galaxy today. 
The first term, which corresponds to the annihilation of particles in a state with L = orbital 
momentum (s-wave), is constrained by helicity arguments to be proportional to (rrif/M x ) 2 , 
and hence very small for light final state fermions. The second term, which corresponds to the 
annihilation in the L = 1 state (p-wave), suffers from the v 2 suppression. For a DM particle 
singlet under the SM gauge group, the radiation of EW gauge bosons from the final state and 
from the internal propagator of the annihilation process eludes the suppressions and opens 
up a potentially sizeable s-wave contribution to the cross section (see Ref. [5] for the case of 
photon radiation and Ref. [6] for gluon radiation). 

In this paper we point out that there is another situation realizing the condition 3 above, 
where we expect therefore the EW corrections to have a great impact. Having in mind the 
Weakly Interacting Massive Particle candidates for DM, it is natural not to restrict oneself to x 
being a gauge singlet, and to consider the possibility that the DM is part of a multiplet charged 
under the EW interactions. In this case, the DM annihilates predominantly in s-wave into 
W + W~, if kinematically allowed. However, now even the initial state of the channel XX ~~ *■ ff 
can radiate a gauge boson; we shall show that this process also lifts the helicity suppression 
and contributes to the s-wave cross section, becoming competitive with the di-boson channel. 



For definiteness, we assume that the DM particle is the electrically-neutral Majorana 
component of a SU(2)l triplet x a > with hypercharge Y = (a wino-like particle) l . The 
coupling with the Z-boson is absent and this DM candidate is compatible with the direct 
detection limits. We neglect the mass splitting of the components of the multiplet, which is 
generated by loop effects [7, 10] and tends to make the charged components slightly heavier 
than the neutral one. The size of this splitting is typically of the order of 100 MeV for a 
TeV-scale DM mass. 

In order to catch the relevance of initial state emission in a model-independent way, we 
work in an effective field theory setup and we restrict to consider the interactions of the triplet 
X a with the SM left-handed doublet L = (f±,f2) T ] the most general dimension-6 operators 
are 

C cS = ^S ab (L 7m P l L) (x V75X b ) + i ^e abc (L ^P L a c L) (x a 7 V) , (1-2) 

where Cd,nd are real coefficients for diagonal and non-diagonal interactions in isospin space, 
Pr,l — n(l i 75) an d & 1 are the Pauli matrices. The assumed Majorana nature of the DM 
forbids some operators that would lead to an s-wave two-body annihilation. The initial state 
radiation lifts the helicity suppression already at the level of dimension-6 operators, unlike 
what happens for the final state radiation, where higher-dimensional operators are needed (see 
Ref. [2] for more details). A more general effective field theory analysis will be presented in 
Ref. [11]. The effective operators in Eq. (1.2) can be generated for instance as the low-energy 
limit of a simple toy model [2, 12] where the DM interacts with the SM left-handed fermions 
through the exchange of a heavy scalar doublet $ 

C mt = -y L La aX a ® + h.c.; (1.3) 

integrating out the scalar sector M$ 3> M x , m w , one obtains the operators in Eq. (1.2) with 
C ND /A 2 = -C D /A 2 = \y L \ 2 /{AMl). 

In the next section we shall discuss the velocity dependence of the amplitudes describing 
the DM annihilation into light fermions, considering both the Final State Radiation (FSR) and 
the Initial State Radiation (ISR) contributions. Subsequently, in section 3 we shall present 
the results for the cross sections and the energy spectra of final particles. Our main results 
are summarized in Section 4, together with prospects for further research. 

2 Velocity dependence of the amplitude 

2.1 Two-body annihilation 

When the DM is part of a multiplet charged under the EW gauge group, and M x > m w , the 
most important two-body annihilation channel is x°X° ~^ W + W~ , which proceeds through 



1 Other representations of the EW gauge group can be considered, see e.g. Ref. [7]. For instance, a simple 
possibility consists of two non-degenerate SU(2)l doublets with opposite hypercharge (higgsino-like particles). 
Direct detection constraints are avoided [7, 8] and also an interesting LHC phenomenology can arise in the 
quasi-degenerate limit [9]. 
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Figure 1: Diagrams for the two-body annihilation into W + W~ {top panel) and the three-body 
annihilation with gauge boson radiation from final {bottom left panel) or initial {bottom right panel) 
state. 

s-wave (see Fig. 1). We shall discuss the importance of this contribution with respect to the 
three-body ISR channel in Section 3 and more thoroughly in Ref. [11]. 

Let us consider instead the annihilation of the DM Majorana fermion into a pair of left- 
handed massless fermions: 

x\k 1 )x\k 2 )^h{pi)Up2). (2.1) 

Only the first term of the effective operators in Eq. (1.2) contributes to this process and 
the annihilation proceeds through the p-wave. The velocity dependence is manifest at the 
amplitude level. In fact, the matrix element for the two-body process is 

1 



M ff~ ■^[u f {p 1 )-f a P L v f {p 2 )} K(fc 2 )7 Q 75U x (fc 1 )] 



(2.2) 



where fc£ 2 = (M x /y/l - {v/2) 2 , 0, 0, ±M x (v/2)/y/l - {v/2) 2 ). The Majorana axial current 
in Eq. (2.2) can be manipulated using the Gordon identities into 



v x {k 2 h a l5U x {h) 



' ' 2 -v x {k 2 h 5 u x {h) - 7 rT r v x (k 2 )(T afi {k 1 - k 2 )p l5 u x {h) . (2.3) 
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For small v, the second term is proportional to v because {k\ — k 2 )^ ~ (0, 0, 0, v M x ); on the 
other hand, in the first term the vector {k\ + k 2 ) a = {p\ +p 2 ) a saturates the final-state current 
in Eq. (2.2) and gives rise to terms proportional to the fermion mass, which are zero in our 
computation. We thus recovered the well-known fact that for Majorana spinors the scattering 
amplitude into massless fermions is proportional to the first power of the relative velocity of 
the incoming particles. 



2.2 Gauge boson emission 

The gauge interactions of the fermion triplet are described by the vector operator e a b c x a W X°-> 
from which it is evident that the neutral component of x does not interact with W 3 . Thus, 



W s is emitted only from the final states, while W can be emitted from either initial or 
final states. Let us study a W^-bremsstrahlung process from either initial or final states, see 
diagrams in Fig. 1. We consider for definiteness the three-body annihilation 

X (h)x°(h) -> fiipi) /aCpa) W~(k). (2.4) 

For the process with FSR, the interaction vertex of two x°'s and two SM fermions is described 
by the isospin-diagonal operator in Eq. (1.2); instead, since ISR changes the isospin state of 
one of the initial legs, the interaction operator is the non-diagonal one. 

Considering the radiation from the final state, the amplitude is the product of the Ma- 
jorana axial-vector current v x 7^75 u x and the fermionic current containing the gauge boson 
emission 
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(2.5) 



where g is the SU{2) L gauge coupling. Using the Gordon identities the Majorana current 
can be simplified as in Eq. (2.3), where the second term is proportional to v, while the 
first term now contains the 4- vector (p\ + p 2 + fc) M ; saturating with the fermionic current in 
(2.5) gives terms proportional to m/, which are zero in our case. Thus, at the lowest level 
in the expansion in M 2 /A 2 , which amounts to restricting to dimension-6 operators, FSR is 
not able to remove the helicity suppression. However, at higher orders (e.g. (9(M 4 /A 4 ) in 
the amplitude) ^-independent terms can arise, for which the inclusion of Virtual Internal 
Bremsstrahlung (VIB) diagrams [5] is crucial. This point is extensively discussed in Ref. [2]. 
On the other hand, the amplitude describing the radiation from the initial state consists 
of the product of the fermionic current Uf j^Pl v / and the Majorana current with the gauge 
boson emission 
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(2.6) 



The effect of ISR is to alter the axial- vector structure of the initial state current, thus pre- 
venting the amplitude from vanishing in the v — > limit 
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(2.7) 



From the amplitudes we studied in this section one can deduce the behaviours of the cross 
sections for the two-body // channel and the three-body channels with ISR, FSR and their 
interference; they can be schematically summarized as 

WffixX ~> fifi) r 
^fsr, isr/fsr(x°X° -> hhW~) r 
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and the precise expressions will be discussed in the next section. The importance of ISR is 
then clear: at the level of dimension-6 operators ISR already opens the s-wave annihilation 
while FSR is still in p-wave. 

3 Results 

3.1 Cross sections 

The cross sections for the various processes are computed from the amplitudes in Eqs. (2.2), 
(2.5), (2.7), with the appropriate coefficients as dictated by the interaction lagrangian (1.2). 
For the two-body annihilation into massless fermions, the cross section reads 

r 2 M 2 

For the three-body processes, it is convenient to define s = (ki + k^) 2 and z = mw/y/s. We 
are working in a situation where the DM mass is larger than the EW scale; so we report 
results as an expansion for small z. Then, the cross sections for the processes with FSR and 
with ISR/FSR interference are 

n 2 C 2 M 2 
vctfsr = f^jf[l5-n 2 + Q\nz(2\nz + 3) + 0(z 2 )]v 2 . (3.2) 

^isr/fsr = - 9 -^§^^[9 + 4lnz + 0(z)]v 2 , (3.3) 

which are both in p-wave. As recalled already in the previous section, an s-wave term in 
the cross section also originates from FSR (together with VIB) but only at a higher order, 
0(M^/A 8 ) in the cross section. Nonetheless, this can have a great impact on the energy 
spectra of final particles [2], but we shall not consider it here. Instead the ISR opens up a 
large s-wave contribution to the total cross section 

9 2 C 2 ND M 2 
w ISR =^^ ir , (3.4) 

which is mediated by the non-diagonal operator in Eq. (1.2) 2 . 

As already anticipated, there is also an important two-body channel x°X° ~~ >" W + W~ , 
whose cross section is not f-suppressed 

WWW = 8^ + ° {zi) ■ (3 ' 5) 

X. 

and can be comparable in size to cxisr. The annihilation channel with ISR can even become 
the dominant one if for instance the W + W~ final state is kinematically forbidden, or if the 



2 The isospin-diagonal operator leads to a spin-dependent elastic DM-nucleon cross section proportional 
to Cq. Experimental bounds from direct detection would therefore constrain Cd, but would not preclude a 
priori the possibility to have a large ctisr- 







Figure 2: Energy spectrum of the emitted W for the three contributions FSR, ISR, ISR/FSR. The 
parameters settings are M x = ITeV, v = 10 -3 . Each distribution is normalized to 1 according to 
Eq. (3.6). 

coefficient Cnd is large. Notice also that the additional g 2 factor in u ww tends to reduce it 
with respect to okr, if Cnd ~ 1; furthermore, the number of colours or families of the final 
state fermions, which we do not consider here, also enhances o"isr with respect to oww- 



3.2 Energy spectra 

Let us now turn to the phenomenological implications of ISR, namely the effect of the W- 
bremsstrahlung on the energy spectra of stable particles resulting from the hadronization and 
decay of the annihilation products. 

In this paper we are interested in extracting the features distinguishing ISR from other 
kind of processes, like the two-body W + W~ channel and the three-body FSR, although their 
cross section can be very different. To make the comparisons more immediate, we define the 
energy spectrum of the channel i as 



dNi 

dE 



1 doi 

Oi dE 



(3.6) 



for i = FSR, ISR, ISR/FSR, WW, so that each contribution is normalized to 1. A more 
detailed analysis will be presented in Ref. [11], where we shall work out an explicit model and 
weigh the various channels appropriately. 

Because the energy of the emitted W boson is entirely distributed among the final particles, 
it is instructive to analyse its energy spectrum separating the contributions from ISR, FSR 
and their interference, see Fig. 2. Notice that the W emission from an ultra-relativistic 
final state particle (FSR) has a characteristic soft/collinear behaviour dN/dy ~ 1/y, where 
y = Ew/M x (see Ref. [1] for further details); on the other hand, the W emission from a 
non-relativistic initial state particle (ISR) shows a somehow peculiar energy spectrum, which 
turns out to be well approximated by the symmetric distribution dN/dy ~ y(l — y). Due 
to the non-relativistic nature of the emitting DM particle and the failure of the factorization 





A" = (E t ,-m f )/M x 



Figure 3: Energy spectra of final positrons (left panel) and antiprotons (right panel), for different 
annihilation channels defined in the text: W + W~ (dashed lines), evW (thick solid lines) and udW 
(thin solid lines). The parameters are set as M x = 1 TeV and v = 10 -3 , and the normalization is 
chosen according to Eq. (3.6). 



property for the three-body cross section, this contribution cannot be caught by the usual 
soft/collinear approximation technique [1]. 

The evolution of the particle species from the primary annihilation products to the final 
stable particles of the SM needs numerical tools. We have carried out this work using our 
own Monte Carlo code, for generating 2—7-2 and 2 — > 3 annihilation events, then interfaced 
to Pythia 8.145 [13] for simulating the subsequent showering, hadronization and decay (see 
Ref. [2] for more details). 

In Fig. 3, we compare the energy spectra of final positrons and antiprotons originating 
from ISR and from WW. We consider two different channels for the ISR case: the lepton 
channel u euW": X°X° ~^ ^L v eiW~ ', eJ / u e iW + ; and the quark channel u udW" : x°X° ~~ >* 
ul(IlW~ , ULdiW + . Interesting features can be extracted from this comparison. For the 
evW channel the very hard positrons can be much more abundant for ISR than for WW, 
due the contribution of the primary positrons. For the udW channel, the antiprotons are 
copiously generated by the W emitted in ISR, especially at low energies because the gauge 
boson is soft, and by the hadronization of the primary quarks, and they can easily overcome 
the antiprotons produced in the WW channel. 

The fluxes of positrons and antiprotons received at Earth can be computed by integrating 
the energy spectra at the interaction point depicted in Fig. 3 over the diffused source consti- 
tuted by the DM distribution in the galactic halo, and then propagating them through the 
halo itself. We do not enter here into any of the details of this process (see e.g. Ref. [14] and 
references therein) but we limit ourselves to a few qualitative considerations. There are irre- 
ducible uncertainties of astrophysical nature that originate from: i) the unknown distribution 
of DM in the halo; ii) the unknown values of the propagation parameters. For the positron 
fluxes at high energy, however, the propagation does not sensibly modify the shape and the 
normalization of the fluxes, for any reasonable choice of the uncertain variables. This is just 
because high energy positrons have anyway originated very close to the location of the Earth 



and therefore have not been affected much during their travel. In turn, this implies that the 
spectral features at high energy introduced by ISR, apparent in Fig. 3, are preserved in the 
final positron fluxes, at least in the case of large DM mass. On the other hand, the fluxes 
of antiprotons are affected by large astrophysical uncertainties both in normalization and in 
shape all across the energy range, making it more difficult to disentangle the different spectral 
shapes. 

4 Summary and outlook 

The inclusion of EW corrections is an essential ingredient to be taken into account for indirect 
DM searches. In this paper we have assumed that the DM is the neutral Majorana component 
of a multiplet charged under the EW interactions and considered the effect of gauge boson ra- 
diation from the initial state of the DM annihilation process. We have restricted ourselves to 
the case where the multiplet containing the DM particle is a SU(2)l triplet, but it is straight- 
forward to work out the cases of multiplets transforming under different representations of 
the EW gauge group. 

The natural annihilation channel for such a candidate is of course through s-wave into 
W + W~ , if kinematically allowed, while the annihilation into light SM fermions is helicity 
suppressed and proceeds through p-wave. However, we found that the iy-bremsstrahlung 
from the initial state removes the suppression and adds a potentially sizeable contribution to 
the s-wave cross section. The gauge boson emission alters the energy spectra of final stable 
particles in a distinguishable way and cannot be ignored for reliable predictions to be used 
for indirect DM searches. 

In a forthcoming paper [11] we shall expand the idea of this work in several directions: a 
complete effective field theory analysis, a full calculation in the context of an explicit model 
which will also allow to weigh precisely the different channels, a detailed computation of fluxes 
of stable particles with the inclusion of the propagation effects. 
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